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Abstract 

lnterleukin-12 is a pro-inflammatory cytokine which promotes Thl and cytotoxic T lymphocyte activities, such as Interferon- 
y secretion. For this reason lnterleukin-1 2 could be a powerful therapeutic agent for cancer treatment. However, Interleukin- 
12 is also excessively toxic, lnterleukin-27 is an immunoregulatory cytokine from the lnterleukin-12 family, but it is not as 
toxic as lnterleukin-12. In recent years, lnterleukin-27 has been considered as a potential anti-tumor agent. Recent 
experiments in vitro and In vivo have shown that cancer cells transfected with IL-27 activate CD8 + T cells to promote the 
secretion of anti-tumor cytokines lnterleukin-1 0, although, at the same time, IL-27 inhibits the secretion of Interferon-y by 
CD8 + T cells. In the present paper we develop a mathematical model based on these experimental results. The model 
involves a dynamic network which includes tumor cells, CD8 + T cells and cytokines lnterleukin-27, lnterleukin-1 0 and 
Interferon-}'. Simulations of the model show how lnterleukin-27 promotes CD8 + T cells to secrete lnterleukin-10 to inhibit 
tumor growth. On the other hand lnterleukin-27 inhibits the secretion of Interferon-}' by CD8 + T cells which somewhat 
diminishes the inhibition of tumor growth. Our numerical results are in qualitative agreement with experimental data. We 
use the model to design protocols of IL-27 injections for the treatment of cancer and find that, for some special types of 
cancer, with a fixed total amount of drug, within a certain range, continuous injection has better efficacy than intermittent 
injections in reducing the tumor load while the treatment is ongoing, although the decrease in tumor load is only 
temporary. 
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Introduction 

Interleukin- 1 2 (IL- 1 2) is a pro-inflammatory cytokine that plays 
a central role in the connection of the innate resistance and 
adaptive immunity, by promoting Thl and cytotoxic T lympho- 
cyte (CTL) activities, such as IFN-y secretion. IL-12 could be a 
powerful therapeutic agent to eradicate tumor or to prevent the 
development of metastasis [1—4]. However, IL-12 has also been 
shown to be excessively toxic [5,6], although there is at least one 
ongoing clinical trial with IL- 1 2 using a new delivery method (IL- 
1 2 DNA plasmid) that is intended to overcome toxicity problems. 
In recent years there has been increasing interest to investigate the 
role of another member of the IL- 1 2 family, namely, Interleukin- 
27 (IL-27), which is less toxic than IL-12, as a potential anti-tumor 
agent [7]. IL-27 is a cytokine capable of regulating Thl, Th2, 
Thl 7, and T re? responses [8]. In autoimmune diseases, Mur- 
ugaiyan et al. [9] have shown that IL-27 promotes production of 
IL- 1 0 and IFN-y by naive human CD4 T cells, and Stumhofer et 
al. [10] demonstrated that IL-27 promotes production of IL-10 by 
CD4 + and CD8 + T cells. Reviewing the role of IL-27 in anti- 
cancer immunotherapy, Swarbrick et al. [11] asserted that IL-27 
may have both pro-inflammatory and anti-inflammatory func- 
tions, and it promotes Thl immune response and CD8 + cell 



activation. Since Hisada et al. [7] first reported the anti-tumor 
efficacy of IL-27 in 2004, the potent anti-tumor activity of IL-27 
has been verified in various tumor models [11-13]. Many studies 
suggest a role of IL-27 in enhancing anti-tumor CD8 + T cell 
responses [7,14-17]. The enhancing role of IL-27 in generating 
anti-tumor CTL response was also demonstrated using IL-27R 
deficient mice [18,19]. 

IL-10 has inhibitory and stimulatory effects on human CD8 + T 
cells [20], and in viral infection it is known to inhibit effector and 
memory CD4 + T cell responses but not memory CD8 + T cells 
[21]. IL-10 may have positive or negative effect on tumor 
suppression (Asadullah et al. [22]). Numerous studies (e.g. [23,24]) 
show that increase in IL-10 produced by macrophages is 
associated with tumor progression, while other studies [25-28] 
suggest that IL-10 plays a positive role in tumor rejection. IL-27 
can induce production of IL-10 in CD8 T cells [10,29]. In a 
recent study, Liu et al. [30] used PI CTL TCR transgenic mouse 
model and mouse plasmacytoma tumor system to investigate how 
IL-27 enhances the anti-tumor responses. They found that IL-27 
significandy enhances the survival of activated tumor antigen 
specific CD8 + T cells in vitro and in vivo, and induces IL-10 
upregulation in these T cells. It was also suggested in [30], and 
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demonstrated in [25-28], that CTL IL-10 production contributes 
to tumor rejection. These results have important implications for 
designing IL-27-based immunotherapy against human cancer. 

In the present paper, we develop a mathematical model that 
describes the anti-tumor activity of CD8 + T cells in terms of IFN-y 
and IL-10 productions, when these T cells are activated by IL-27 
from the tumor microenvironment. The model is based on the 
experiments by Liu et al. [30] (with mice infected with 
plasmacytoma) whereby cancer cells are transfected with an IL- 
27 vector so that IL-27 is released in the tumor microenvironment. 
We show that the model simulations agree qualitatively with the 
experimental results of [30] . We next extend the model to include 
therapeutic treatment of cancer in wild type mice by injection of 
IL-27. We note however, that in this case, the tumor microen- 
vironment includes both CD4 + and CD8 + T cells (whereas in the 
experiments with transgenic mice of [30] the CD8 + T cells were 
taken to be dominant). As mentioned above, IL-27 promotes the 
secretion of IL-10 and IFN-y by CD4 + T cells [9-11], and we 
assume that these cytokines have the same effect on tumor 
rejection as those secreted by CD8 + T cells. We then use CD8 + T 
cells to represent both cells, CD4 + and CD8 + . In the modified 
model the only source of IL-27 comes from the drug, since cancer 
cells do not generally secrete IL-27. We compare the efficacy of 
different strategies of IL-27 injections. For example, we found that 
continuous injection of IL-27 for 24 weeks at a fixed amount F, 
within a certain range, is more effective than intermittent injection 
of the amount IF , full three weeks at a time with three weeks 
spacing between injections, for 24 weeks. These predictions 
however must be viewed just as suggestions since they may only 
apply to special types of cancer, such as plasmacytoma in bone or 
soft tissue, and since, furthermore, the model does not include 
other important factors in tumor progression such as angiogenesis 
and the immune response. 

Results 

Mathematical model 

In this model, we assume that the tumor is spherical and that it 
initially lies in a spherical tissue of radius L. The variables that will 
be used in the model are listed below and we assume that all the 
variables are radially symmetric: 



Tumor cell (c) 
secrete 



h 



IL-10 (7 I0 ) 




t IL-27 (7 27 ) N 


activate 





' immune 
response 



c^Tcell Cr)-j5555 |gN-r(/,,) 



Figure 1. A network of IL-27. A network showing how IL-27 affects 
the immune response to tumor cells. CD8 + T cells are activated by P1 A 
antigen from tumor cells as well as by IL-27 which is secreted by tumor 
cells. Activated CD8 + T cells secrete IFN-y, inhibited by IL-27, and IL-10 
enhanced by IL-27. IL-10 and IFN-y inhibit tumor cells. 
doi:1 0.1 371/joumal.pone.0091 844.g001 

partial differential equations based on the network exhibited in Figure 1 . 
The parameter values are estimated in Methods. In our model we shall 
include diffusion of cells and cytokines, as was done in many other 
models of solid tumors (which include plasmacytoma [31-38]) 



IL-27. The following equation describes the evolution of 
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The first term represents the diffusion of hi with coefficient 
D/, 7 . Although diffusivities of cytokines and cells may depend on 
the concentrations of the tumor cells and normal healthy cells, for 
simplicity, here and in the sequel all diffusivities are assumed to be 
spatially uniform. In the experiment in [30], Liu et al. used gene 
transfected tumor cells, J558-IL-27, to produce hi m the tumor 
microenvironment. Accordingly, we use the second term to 
describe the production of 727 by the transfected J558-IL-27 
tumor cells. The last term stands for the degradation oi hi- The 
parameter values of Equation (1) are given in Table 1. 

IL-10. The Interleukin- 1 0 (IL-10) in Figure 1 is pro-inflam- 
matory, in accordance with the experiments of [30] . It satisfies the 
equation: 



hi(r,t) 
ho(r,t) 
T(r,t) 

Iy(r,t) 
c{r,t) 



concentration of Interleukin — 27 in pgjcm 3 ' , 
concentration of Interleukin — 10 in pg/cm 3 , 
(tumor antigen specific) activated CD8 + T 
cell density, cell /cm 3 , 

concentration of Interferon — )' in pg/cm 3 , 
tumor cell density, cell /cm 3 , 



where r is the distance from a point X = {X\,X2,xi) to the origin: 
r = \ x\ + Xj + x^ and 0 < r < L. These variables satisfy a system of 



5/ io n 1 8 , 37io 
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production by CTL promoted by /27 
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Table 1. 


Parameters for the IL-27 equation. 








Parameter 


Description 


Value with unit 


Reference 




diffusion coefficient of I 2 j 


1.25 x l0- } cm 2 /day 


[48] & estimated 




production rate of I 21 from tumor 


1.5 x 10- 5 pg/ cell /day 


[30] & estimated 




degradation rate of In 


2/ day 


[48] & estimated 


doi:1 0.1 371 /journal.pone.0091 844.t001 
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Table 2. Parameters for the IL-10 equation. 




Parameter 


Description 


Value with unit 


Reference 




diffusion coefficient of /io 


1.25 x 10- 3 cm?/day 


[48] 


■Ml) 


production rate from CTL without IL-27 


%.$9xl0- s pg/cell/day 


[30] & estimated 


aio 


max production rate from CTL with IL-27 


I.128xl0- 4 pg/cell/day 


[30] & estimated 


n id 




5 x Mftpg/em* 


[49] & estimated 




degradation rate of I iQ 


l.6xW/day 


[48] 


doi:1 0.1 371 /journal.pone.0091 844.t002 



The first term is the diffusion of I\o. The second term accounts 
for the production of Iio by CD8 + T cells for the absence of ^27 
[30]. The experiments in [30] indicate that /j7 significantly 
increases the production of Iio by CD8 + T cells, and this is 
accounted by the third term. The last term is the degradation of 
Iio. The parameter values of Equation (2) are listed in Table 2. 

CD8 + T cells. The equation for the density of (activated) 
CD8 + T cells, T(x,t), is given by 



dT 



D 1 8 2 8T 

r 2 dr dr 



St 



c+c t 



l+a T I 21 +^ T a T Ii 



immune response 



T. 



(3) 



pro— survival by Ij-] and I^q 

The first term is a dispersion of CD8 T cells with coefficient 
Dj. The second term accounts for activation of CD8 + T cells by 
P1A antigen from the tumor cells. In promotes survival of CD8 + 
T cells, and so does also /in, but to a smaller degree [30]. We 
present these two facts by correspondingly decreasing the death 
rate fi T of T cells in the last term of Equation (3). The parameter 
values in Equation (3) are given in Table 3. Although I%j 
contributes more than Tin to promote the half-life of CD8 + T cells, 
we take /? r = 9 since the concentration of I\o is much smaller than 
the concentration of 727- 

IFN-y. Interferon-y (IFN-y) is a cytokine with diffusion 
coefficient D y and degradation rate It is produced by T cells 
and, as shown in [30], the production is inhibited by 727- Thus, 
I y (x,t) satisfies the equation: 



8I y 
~8t 



D 



1 8 , , 8L 



(r-s L ) + 



a v T- 



production by CTL inhibited by ^7 



' r 2 dr dr 

diffusion 
degradation 



Table 4 lists all parameter values of (4). 

Tumor cells. The density of tumor cells, c(x,t), satisfies the 
following equation: 



8c 



D, 



1 8 , ,3c. 



(rV) + ^(l--) 



r 2 dr dr 



diffusion 

ho 

>lc j T~I c 




(5) 



inhibition by / 10 inhibition by IFN— ) 



The second and third terms represent the proliferation and 
death of cells, respectively. Generally, I\q is regarded as an anti- 
inflammatory cytokine. However, in different experimental 
models, Zio could suppress or promote the functions of immune 
system [30,39] . Liu et al. [30] found that the Zio produced by CTL 
contributes to tumor rejection. Hence, the fourth term accounts 
for the indirect inhibition of tumor cells by /in. Cytokine, T y 
promotes the anti-tumor response, such as increase production of 
IL-12, and induces natural killer cells to kill cancer cells [40,41]. 
For simplicity, we take the fifth term in (5) to represent the 



Table 3. Parameters for CD8 + T cell equation. 




Parameter 


Description 


Value with unit 


Reference 


D T 


diffusion coefficient of CTL 


4.32 x l0- 6 cm 2 /day 


[48] 


St 


production rate of CTL activated by tumor 


1.3968 x 10 s cell / cm 3 /day 


[30] & estimated 


Ct 




5.76 x 10 1( W//W 


estimated 


(7 j 




2 x 10~ i pg/cm > 


estimated 


fir 




9 


estimated 


f'T 


death rate of CTL 


3x 10-' /day 


[48] 


doi:1 0.1 371 /journal.pone.0091 844.t003 
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Table 4. Parameters for the IFN-y equation. 





Parameter 


Description 


Value with unit 


Reference 


Dy 


diffusion coefficient of I y 


1.25 x \f)- y cm 2 /day 


[48] & estimated 


Xy 


max production rate of I r from CTL 


3J2x\0- 5 pg/cell/day 


[30] & estimated 


Sy 




5 x ICfipg/cni 3 


estimated 


t'y 


degradation rate of 


Lib/day 


[49] 



doi:1 0.1 371 /joumal.pone.0091 844.t004 



(indirect) inhibition of tumor cells by Iy. The parameter values are 
listed in Table 5. 

The dimensional and dimensionless values of all the parameters 
of Tables 1-5 are listed in Table 6. 

Initial conditions. We assume that tumor cells are initially 
concentrated near r = 0, taking 

e(r,0) = CQ(e~ r/t — e~ L/e ), where L=l cm or L=1.5 cm, 

and £ is a positive number less than or equal to 1. In the 
simulations, we shall take e = 1 but the results do not change 
qualitatively with smaller values of e. Since Ixi is produced by 
J558TL-27 tumor cells, the initial concentration of hi should be 
similar to the density of tumor cells; we take 

hi(rfi)=—c,(e- rli - e - L l% 
Initially, there are no activated CD8 + T cells, hence 
T(r,0) = 0. 

Since Iio and I y are produced by CD8 + T cells, we take 
7 10 (r,0)=0 and I r (r,0) = 0. 

Boundary conditions. Since all variables are radially sym- 
metric, the first r-derivative at r = 0 is equal to zero. We assume 
no-flux condition for all variables at r = L. This is justified by the 
fact that L is large enough so that the exterior of the ball of radius 
L lies completely within the healthy tissue, initially. T = 3 days, L 0 = 0.5 cm 



Table 5. Parameters for tumor cell equation. 




Parameter 


Description 


Value with unit 


Reference 




diffusion coefficient of tumor 


8.64 x W- 6 cm 1 /day 


[48] 


h 


max proliferation rate 


4.68 x 1Q- 1 /day 


[30] & estimated 


ft. 




10 s cell/cm 3 


[48] 


V-c 


death rate of tumor 


1.73x10-' /day 


[48] & estimated 


1c 


inhibition rate of tumor from 7i 0 


3.45 x 10-' /day 


[30] & estimated 


Oc 




1.5 x Xtfpg/cm 1 


estimated 


If 


inhibition rate of tumor from IFN-y 


6x 10-' /day 


[30] & estimated 


Sc 




3 x I0 1 pg/cm s 


estimated 


doi:1 0.1 371 /joumal.pone.0091 844.t005 
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Parameters nondimensionalization. We nondimensiona- 
lizate the Equations (1) - (5): 

f = r/Lo,i = t/r, 

h7 = h7/l2 1 M 0 = I W /Ii 0 ,f=T/T°J r =L/I°;c = c/c 0 , 
^0 

{a 27 Ao.Oy} = t:{co«27 I 'I$ 7 ,T 0 aiio/ 'I° 0 ,T 0 x.,/ 7"} , 
{Az7.Aio.Ar.Aj. A} = ^{fe7.Mio./tr./^,Mc}. 
{s m ,s T } = t{ r 0 Jio/iio,ffr/ 7o},{^,.,%} = z{r\ c ,riy}, 

{&io,S y } = 4j- {a w ,Sy},a T = 4i°t,Pt = Pt, 

a c = a c /I° 0 ,s c = s c /I°,2.i = xl\ ,c t = c t /c 0 ,c T = c T / cq, 
where 
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Table 6. Model variables and units. 



Parameter 


Dimension value 


Dimensionless value 


D hl 


1.25x10- 3 cm 2 /day 


1.5 x 10" 2 


Di,, 


1.25 xlO" 3 cm 2 /day 


1.5x 10~ 2 


D T 


4.32 x lO-' 1 cm 1 / day 


5.184x 10- 5 


D, 


1 .25 xlO- 3 'cm 2 /day 


1.5x 10- 2 


D e 


8.64 x 10 -6 'cm 2 /day 


1.0368 x 10~ 4 


f'n 


2 / day 


6 


f'w 


1.6 x 10/ 'day 


4.8 x 10 


f'T 


3x10"' /day 


9x 10-' 


Hy 


2.16 /day 


6.48 


f'c 


1.73x10-' /day 


5.19x10-' 


a-n 


1. 5 x 10- 5 pg/ cell /day for 
small production 


3.24 x 10 2 




7.5 x 10~ 5 pg/ cell /day for 
moderate production 


3.24xl0 2 x5 




3 x 10" A pg/ cell /day for large 
production 


3.24xl0 2 x20 


y.w 


1.128 x 10-" pg/ cell /day 


3.384x 10-' 


0-7 


3.72xl0- 5 pg/cell/day 


1.116x10-' 


Pt 


9 


9 


Sw 


8.89 x 10- % pg/ cell /day 


2.667 x 10-" 


St 


1.3968 x 10 8 cell/ cm 3 / day 


4.1904 x 10 3 


s r 


5 x 10 G pg/cm 3 


5x 10 4 


$c 


3 x 10 2 pg/cm 3 


3 




5 x lO^pg/cm 


5x 10 


ft J 


2xl0- 4 pg/cm 3 


2x 10- 2 


a c 


1.5 x 10 2 pg/cm 3 


1.5 


1c 


3.45 x 10-' /day 


1.035 


>h 


6x 10-' /day 


1.8 


h 


4.68 x 10-' /day 


1.404 


c. 


Wcell/cm 3 


1 .38889 


Ct 


5.76 x 10">cell/cm } 


8x10 


L 


1cm in Figs. 2-5 


2 




1.5cm in Figs. 6-10 


3 


X 


Mays 


1 


L 0 


5 x lO-'cm 




CO 


7.2 x 10 s cell/cm 3 


1 


1 0 


i u ceil I L in 






lO^-pg/cn? 




To 


Wpg/cm 3 


^^^^^^^^ 


f 


10 2 pg/cm 3 




F 


3 x lOftpg/cm 3 /day in 
Figs. 6, 7, and 10 


9x 10 2 




1 x lO^pg/cm 3 /day in 
Figs. 8A and 9A 


3x 10 2 




0.5 x 10 4 pg/cm 3 /day in 
Figs. 8B and 9B 


1.5 x 10 2 


a 


2.25cm 2 


9 
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c 0 = 7.2 x 10 s cell /cm 3 , T 0 = \0 5 cell/c 



4 1 =P l0 =Py=l&pg/cm\ 



For nondimensional variables and parameters, we consider the 
tumor growth in a ball {f<2} or {r<3}. The nondimensional 
PDE model is given by the following system of equations: 



5727 

di 



ei 



dT 

77 



dly 

7F 



dc 
di 



27 f2 gf flf 



027C - £27/27 

production by tumor degradation 



bj — — (f 2 — ) + 

10 j.2 flj. flf 



diffusion 

+ aioT- 

production by CTL promoted by ly-* 

1 8 ,„ 2 df x . „ c 




r L or or 



c + cr 



diffusion 



1 + arhi + PtQtIio 

v ' 

pro— survival by /07 anc < 7|q 

1 8 „. dly 



immune response 

t 



(6) 



v f 1 df ^" 8r ^ 



(XyT- 



production by CTL inhibited by lyn 



degradation 

z>4^ 2 ?) + ^ ( i-A) 

r z or or c* 



proliferation 

>h~ -c 



death 



- ho . 

*7 C ~ c 

IlQ + dc " Iy+Sc 

inhibition by I\q inhibition by IFN— y 



Numerical simulation 

The model (6) was simulated, in nondimensional variables, 
using madab with dr = 0.02 and di = 0.0l (i.e., dr = 0.01cm and 
dt= j^day in dimensional units). Four cases were considered: 

(i) J558-Ctrl tumor cells. 

(ii) J558TL-27 tumor cells with small production rate of IL-27. 

(iii) J558-IL-27 tumor cells with moderate production rate of IL-27. 

(iv) J558-IL-27 tumor cells with large production rate of IL-27. 
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It has been reported in [30] that In can enhance the population 
of CD8 + T cells. Moreover, 727 also enhances 7io produced by 
CD8 + T cells to inhibit the tumor growth, but at the same time it 
suppresses the pro-inflammatory cytokine I y secreted by CD8 + T 
cells. In spite of its inhibition of 727 still promotes CD8 + T cells 
to suppress the tumor growth. 

In view of these experimental results we expect the total mass of 
7io to increase from cases (i) to (iv), the total CD8 T cell 
population to increase from cases (i) to (iv), and the total 
population of cancer cells to decrease from cases (i) to (iv), as 
time progresses. 

Correspondingly, we associate with the four cases (i) - (iv) 
increasing values of the parameter GC27: 

a 2 7 = 0 and 7 2 7(r,0) = 0,re[0, 2] in Figures. 2 — 5, 
in case (i) of J558 — Ctrl tumor cells, 

<x 2 7 = 0 and ln(r,0) = 0,re[0, 3] in Figures. 6 — 9, 
in case (i) of J558 — Ctrl tumor cells, 

a27 = 324(i.e., (X27 = 1.5 x 10~ 5 pg/cell/day), in case (ii), 

6t 2 7 = 324x 5(i.e., a. 2 i = 1.5x\0~ 5 pg/cell/day), in case (iii), 

6t27 = 324 x 20(i.e., OC27 = 3 x pg/cell/day), in case (iv). 



In Figures 2-5, we took L= \ cm such that tumor cells are not 
visible near the boundary r=\cm for all time te[Q,90days]. 
Figure 2 shows the time-dependent profiles of the total mass of 



InJiaJy, and total populations of CD8 + T cells and cancer cells. 
We see that growth/decrease of these variables, as 0.21 varies, 
corresponds qualitatively to the experiments in [30]. Figures 3-5 
show significant spatial variations of these variables at days 3,9, 
and 15, with or without production of 727. We also see the effect of 
^27 on cancer cells and CD8 + T cells densities at different distances 
from the point of origin of the cancer. For example, at the origin, 
at day 3 the cancer cells density changed from 4.95 x Xif 1 cell / cm 3 
with no treatment by ^27 to 4.65 x cell / cm 3 with largest 
production of ^27, while at day 15 it changed from 
4.2 x 10 s cell /cm 3 with no treatment to 1.7 x W^cell/cm 3 with 
largest production of 727- Similarly, at the origin, the CD8 + T cell 
density increased at day 3 from 2.3 x \0 6 cell /cm 3 without 
treatment to 3.3 x 10 6 cell /cm 3 with largest production of 727, 
whereas at day 15 the density increased even more significandy 
from 3.5 x W 6 cell/cm 3 with no 7 2 7 treatment to 9 x W 6 cell/cm 3 
with maximal production of 727- Note that Figure 5E and 5F show 
that the tumor cell density is almost zero near the boundary 
r= 1cm and the tumor cells concentrate in the region r = 0.5cm. 

Tumor initiating in internal organs can also be treated by 727, 
but the mechanism for introducing 727 will depend on the location 
of the tumor. For example, in colitis induced colon cancer, one 
could use yeast which were programmed to express 727 [4-2]. 
Oncolytic virus which are engineered to produce 727 within tumor 
cells could turn the tumor into immunogenetic, thus enabling the 
immune system to reject the tumor. 

We want to use our model in order to design treatments for a 
wild type mouse by 727 injection. We recall, as noted in the 
Introduction, that for wild type mouse, both CD4 and CD8 + T 
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Figure 2. Evolution of cells and cytokines for different production rates of IL-27. (A), (B), (C), (D), and (E) are the profiles of total number of 
hi,ha,T,I y , and c, respectively, within 90 days. In (E), the curves displayed from top to bottom are for J558-Ctrl tumor cells, J558-IL-27 tumor cells 
with small («27 = 1.5 x W~ 5 pg/ cell /day), moderate (0(27 = 7.5 x \Q~ 5 pg/ cell /day), and large [_azi = '5xl0~*pg/cell/day) production of IL-27, 
successively; L= 1cm. 
doi:1 0.1 371 /journal.pone.0091 844.g002 
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Figure 3. Spatial distributions at day 3. (A), (B), (C), (D), and (E) are the spatial distributions of Ixi,ho,T,I y , and c, respectively, at day 3 for 
different production rates of IL-27. (F) is zoomed in (E) near r=0.5cm. In (E), the curves displayed from top to bottom are for J558-Ctrl tumor cells, 
J558-IL-27 tumor cells with small (227 = 1.5 x pg/cell/day), moderate (3.27 = 7.5xl0~ 5 pg/cell/day), and large (a.xi^xX^ 4 pg/cell/day) 
production of IL-27, successively; L=\cm. 
doi:1 0.1 371 /journal. pone.0091844.g003 



cells produce IL-10 and IFN-y [9-11] and we assume that IL-10 
secreted by CD4 + T cells has the same tumor rejection quality as 
the IL-10 secreted by CD8 + T cells. We then use CD8 + T ceUs to 
represent both cells, CD4 + and CD8 + . We also note that in vivo 
tumor cells do not generally secrete ^27, so we take 0:27 = 0 in 
Equation (1). But we also need to include an injection term in 
Equation (1) for Ixi. If we denote the injection density by f(r,t) 
then Equation (1) becomes 



2| for 0<r<L and 0<?<24 weeks 

. L +" (in real time) , 

' { ' J) ~ ' for 0<r<L and 24</<30 weeks ^ ' 

0, 

(in real time) , 
and for the intermittent strategy 



3/27 „ 10,2 5/ 2? ^ , /y \ 



- Ihnhl ■ ( 7 ) 

injection of ^7 degradation 



We make the pharmacokinetic assumption that f(r,t) decreases 
in r from the outer boundary of the tumor (r= 1.5 cm) towards the 
inner core (r = 0), and take 



f(r,t) = Fx 



L 2 + c 



(8) 



where a is some positive constant; F is viewed as the "amount" of 
injection. 

We consider here, for illustration, two strategies of treatment: (i) 
continuous injection of 727 a t a fixed amount F for 24 weeks, and 
(ii) intermittent injections, at double amount IF, full three weeks at 
a time with three weeks spacing between injections. Accordingly, 
for the continuous strategy 



2Fx 4±«, 

iA+a 



0, 



for 0<r<L and ?2i<?<fe+i 
i =0,1,2,3 

for 0<r<L and ? 2 / + i < t< ?2(/+i), v 
i= 0,1,2,3 and 24</<30 weeks 



(10) 



in case (ii), where the length of each interval is three weeks 

(the drug is injected only during the intermittent intervals 
teij'a+l]) an d ?o = 0. In the following simulations, we take 
fl = 2.25cm 2 ; however the same results remain qualitatively the 
same for other values of a (not shown here). 

In Figures 6-9, we take L = 1. 5cm so that the tumor cell density 
remains negligible near the boundary r = L, during the entire 
simulation time which is 30 weeks and hence the boundary 
conditions are not affecting the results during the entire simulation 
(For longer simulation time, e.g. 60 weeks, we need to take 
L = 3cm (not shown here.)). We also take the simulation mesh size 
dr = Q.Q\cm and dt = 0.0lday. In Figure 6, we compare the results 
of the two strategies in case F = 3 x \0 4 pg/cm 3 /day. We see that 
continuous injection yields better results in reducing the tumor 
level and slightly delaying relapse after the drug is withdrawn. 
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Figure 4. Spatial distributions at day 9. (A), (B), (C), (D), and (E) are the spatial distributions of I 2 7,I\o,T,I r , and c, respectively, at day 9 for 
different production rates of IL-27. (F) is zoomed in (E) near r=0.5cm. In (E), the curves displayed from top to bottom are for J558-Ctrl tumor cells, 
J558-IL-27 tumor cells with small (227 = 1.5 x W~ 5 pg/cell/day), moderate (3.2i = 7-5xl0~ 5 pg/cell/day), and large (3.27 = 3xl0~ 4 pg/cell/day) 
production of IL-27, successively; L=\cm. 
doi:10.1371/journal.pone.0091844.g004 
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Figure 5. Spatial distributions at day 15. (A), (B), (C), (D), and (E) are the spatial distributions of In,I\o,T,I y , and c, respectively, at day 15 for 
different production rates of IL-27. (F) is zoomed in (E) near r = 0.5cm. In (E), the curves displayed from top to bottom are for J558-Ctrl tumor cells, 
J558-IL-27 tumor cells with small (227 = 1.5 x W~ 5 pg/cell/day), moderate (3.2i = 7-5xl0~ 5 pg/cell/day), and large (3.27 = 3xl0~ 4 pg/cell/day) 
production of IL-27, successively; L=\cm. 
doi:10.1371/journal.pone.0091844.g005 
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Figure 6. Comparison of continuous versus intermittent treatment. (A), (B), (C), (D), and (E) are the profiles of total number of hi,ho,TJ y , 

and c, respectively, for model (6) with L= 1.5cm which the first equation for /?7 is replaced by (7) and all parameter values are taken from Table 6. In 
(E), the upper curve is for J558-Ctrl tumor cells, the dotted-dashed curve (—• — •—) is for intermittent injection, and the dashed curve is for 
continuous injection with F = 3 x 10 4 pg j cm 7, / day and a = 2.25cm 2 , for the first 24 weeks. 
doi:10.1371/journal.pone.0091844.g006 
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Figure 7. Concentration profiles of tumor cells at different times. (A), (B), and (C) are the concentration profiles of c at times 4 weeks (short 
time), 24 weeks (time at which injections are withdrew), and 30 weeks (the final time for simulation), respectively, under drug amount 
F = 3x 10 4 pg / en? / 'day and a = 2.25cm 2 . The upper curve is for J558-Ctrl tumor cells, the dotted-dashed curve ( — ■ — ■—) is for intermittent injection, 
and the dashed curve is for continuous injection. The concentration of tumor cells are not visible, when r is close to 1.5cm, for all te{0,30weeks]; 
L= 1.5cm. 

doi:1 0.1 371 /journal.pone.0091 844.g007 
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Figure 8. Comparison of continuous versus intermittent treatment for different drug amount. (A) and (B) are the profiles of total number 
of c with F — 1 x W 4 pg/cm 3 /day and ^ = 0.5 x /day, respectively, for model (6) with L= 1.5cm and a = 2. 25cm 2 which the first equation 

for hi is replaced by (7) and all parameter values are taken from Table 6. The upper curve is for J558-QH tumor cells, the dotted-dashed curve 
(—• — ■—) is for intermittent injection, and the dashed curve is for continuous injection. 
doi:10.1371/journal.pone.0091844.g008 
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Figure 9. Tumor growth and migration for shorter injection schedule. (A) is the concentration profile of c at 4 weeks and (B) is the profile of 
total number of c for te[0,4weeks], under f = 3 x 10 4 pg/cm i /day, for model (6) with L= 1.5cm and a = 2.25cm 2 which the first equation for /27 is 
replaced by (7) and all parameter values are taken from Table 6. The upper curve is for J558-Ctrl tumor cells, the dotted-dashed curve (—■ — •—) is for 
intermittent injection, and the dashed curve is for continuous injection. (A) shows that the concentration of tumor cells are not visible near the 
boundary r= 1.5cm, for all te[0,4weeks]. 
doi:10.1371/journal.pone.0091844.g009 
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Figure 7, for the same experiment as in Figure 6, shows the 
concentration profiles of tumor cells at times 4 weeks, 24 weeks, 
and 30 weeks for J558-Ctrl, intermittent injection, and continuous 
injection cases. Notice from Figures 7A, 7B, and 7C that the tumor 
has progressed during the periods of 4 weeks, 24 weeks, and 30 
weeks to r = 0.6cm, r= 1.05cm, and r= 1. 2cm, respectively. 
Figures 6 and 7 show that hi injection slows down tumor growth 
during drug injection, but it does not change the migration speed 
of tumor cells. Figure 8 compares the results of the above two 
strategies for smaller values of F, namely, F=lx 10'' pg / cm / day 
and F = 0.5 x I0 4 pg/cm 3 /day. We see that continuous injection is 
still more effective, but, for smaller amount of injection, the 
relative advantage of continuous injection is decreased. Simula- 
tions of these two strategies for other values of F in the range of 
10 3 — \0 6 pg/cm 3 /day (not shown here) give the same results, 
namely, that continuous injection is preferable to intermittent 
injections. In order to make a definite recommendation on 
continuous versus intermittent injection one would need to 
consider also possible side-effects that may arise from these two 
strategies. 

Although the expected lifespan of the mouse in the experiments 
of Liu et al. [30] was one month, for the purpose of therapy we 
performed simulations for the longer period of 30 weeks. But it is 
also interesting to consider the case of treatment for one month 
only. This is done in Figure 9 where we have taken in Equation (9) 
0 < t < 4 weeks for the continuous treatment, and, in Equation 
(10), intermittent time tj=j weeks. Figure 9 shows that continuous 
treatment is again preferable to intermittent treatment. Figure 9A 
shows the concentration profile of tumor cells at 4 weeks; note that 
tumor cells do not reach the boundary within 4 weeks. Figure 9B 
displays the profile of total number of tumor cells. We see that the 
continuous treatment still has better efficacy than intermittent 
treatment. 

Sensitivity analysis 

In order to provide support to the robustness of the simulation 
results, we ran sensitivity analysis on parameters which appear in 
Equations (1) - (5). The parameters chosen for the sensitivity 
analysis are either those whose baseline was crudely estimated, or 
those that seem to play more important role in the model 
predictions. 

We list these parameters with their ranges, baselines, and units, 
in Table 7. In this analysis, 0:27 varies from 32.4 to 12960. 
Following the sensitivity analysis method described in [43], we 
performed Latin hypercube sampling and generated 5000 samples 
to calculate the partial rank correlation coefficients (PRCC) and p- 

value, with respect to the ratio J„ 2 c(r,t)r 2 dr/ Jj 2 c~ (r,t)r 2 dr, where 
c(r,t) (resp. c~(r,t)) accounts for theJ558-IL27 (resp. J558-Ctrl) 
tumor cell density, at ? = 30 and re[0,2]. The PRCC and their p- 
values are listed in Table 8. A negative PRCC (i.e. negative 
correlation) means increase in the parameter value will decrease 
the ratio Jg c(r,i)r 2 dr/ J 2 c~ (f,t)f 2 dr; that is, it will increase the 
rejection of tumor treated by IL-27 versus untreated tumor. 
Conversely, positive PRCC means that increased rejection of the 
tumor (treated by IL-27 versus untreated) will occur if this 
parameter is decreased. 

The sensitivity analysis data are shown in Figures S1-S4 in 
Supplementary Material and summarized in Table 8. The most 
significant negatively correlated parameters in promoting rejection 
of tumor treated by IL-27 versus untreated tumor are Hxi^t^t', 
less significant parameters are r\ y ,r] c ,a.\$,fi T . The effect of (X27 has 
already been displayed in Figures 2-5. The negative correlations 
of St, '7 c , and 17,, are not surprising, since St is the rate by which 



Table 7. Parameters chosen for sensitivity analysis. 





Parameter Range 


Baseline 


Unit 



«27 

am 
st 
■5" 1 0 

n c 
% 

s ( 

Sy 

h 

C T 

doi:1 0.1 371 /journal.pone.0091 844.t007 

tumor activates T cells (while T cells are increased with IL-27 
treatment; see Figures 2-5) and v\ c and f\„ are, respectively, the 
killing rates of tumor cells by I\o and (while I\o and I y increase 
with IL-27 treatment; see Figures 2-5). The negative correlations 
of aio,(7r, and fi T are also not surprising, since ^27 promotes the 
production of I\a to inhibit tumor cells, and larger oj and f> T 
promote survival of CD8 + T cells. 

The most significant parameters in promoting tumor are 
g c ,Ct,s c and, to a smaller degree, G\q. This also is not surprising, 
since increasing a c results in decreased inhibition of tumor cells by 
iio, increasing cj results in decreased number of CD8 + T cells, 
increasing s c results in decreased inhibition of tumor cells by I y , 
and increasing rj 10 results in decreased Iio- We note that the 
parameters Sin and Sy, in Table 8, have small PRCC with p-values 
that are larger than 0.01; this means that they are not sensitive to 
the ratio J Q 2 c(r,i)r 2 dr/ J 2 t" (f \t)r 2 df . 

Table 8. The PRCC and p-value of parameters for sensitivity 
analysis. 



Parameter 


PRCC 


p-value 


a 27 


-0.63709 


<0.01 


S10 


0.0074315 


0.59933 


S10 
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-0.52728 
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n c 


-0.30727 
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-0.38691 
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Sy 


0.010033 


0.47817 


(Tf 


-0.59115 


<0.01 


Pt 


-0.10923 


<0.01 


Ct 


0.52848 


<0.01 
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Discussion 

IL-12 plays a central role in linking the innate resistance and 
adaptive immunity, and could be a powerful anti-tumor agent. 
However, since IL-12 is excessively toxic, the cytokine IL-27, 
which is a less toxic member of the IL-12 family, has been 
considered as a possible replacement of IL- 1 2 as anti-tumor agent 
[7,8,14-17]. It was demonstrated by Liu et al. [30] that IL-27 
enhances the survival of tumor antigen specific CD8 + T cells and 
induces their upregulation of IL- 1 0, which acts as an anti-tumor 
cytokine. This suggests that IL-27 could play an important role in 
immunotherapy against human cancer. 

The aim of the present paper was to develop a mathematical 
model that can be used to explore and predict the efficacy of 
different protocols of IL-27 treatment. To do that we first set up a 
dynamical system of partial differential equations whereby IL-27 is 
produced by transfected J558-IL-27 tumor cells, as demonstrated 
in the experiments of Liu et al. [30]. The model included IL-27- 
induced CD8 + T cells and cytokines IL- 1 0 and IFN-y. By carefully 
estimating the parameters of the equations we showed that the 
model simulations agree with the experimental data of Liu et al. 
[30]. 

The model can be used to examine the effect of injecting IL-27 
into the microenvironment of cancer in a mouse, and design 
strategies for such injections. We illustrated this by comparing the 
efficacy of two protocols: (i) continuous injection (e.g., daily) of IL- 
27 for 24 weeks at a fixed amount F, and (ii) intermittent injections 
during the first 24 weeks with three weeks injection at a fixed 
amount 2F followed by three weeks spacing, and withdrawing the 
drug after the 24 weeks for both protocols (i) and (ii). We found 
that the continuous injection has better efficacy in reducing the 
tumor load, and also in delaying relapse after the drug is 
withdrawn, while the treatment is ongoing. However, in estab- 
lishing these results we made the assumption that IL-10 produced 
by IL-27 activated CD4 T cells has the same pro-inflammatory 
property as the IL-10 produced by CD8 + T cells. In addition, we 
made the pharmacokinetic assumption that the drug density 
decreases toward the inner core of the tumor, and we also took the 
drug "amount" F in the range of 10 3 -W 6 pg/cm 3 /day. 

We note that our model was based on the experiments by Liu 
et al. [30] with plasmacytoma, but not with other tumor cells. 
Furthermore, the model did not include the effects of lymphoid 
and vascular compartments, as these were not reported in [30]. 
Hence the present paper should be considered only as an initial 
building block for a more comprehensive model which should 
include angiogenesis as well as the immune response of 
macrophages, dendritic cells, and T cells (Thl, Th2, Thl7, and 
T r< , g s). We note in particular that pro-inflammatory macrophages 
(Mj) secrete a family of IL-12 cytokines including IL-27 [44], and 
the IL-12 family attracts CTLs which kill tumor cells, so that M\ 
macrophages suppress tumor growth. On the other hand, anti- 
inflammatory macrophages (A/2) secrete IL-10 which promotes 
tumor growth [23,24]. Regulatory T cells promote tumor growth 
and are inhibited by IL-27 [45,46]. Thus the present paper's 
prediction of the efficacy of different protocols of treatment of 
plasmacytoma in bone or soft tissue with IL-27 will need to be re- 
examined when more data become available that will enable us to 
include the important compartments of the immune and vascular 
systems. 

We note also that the proposed intervention with IL-27 in our 
paper shows benefits only while the treatment is ongoing. The 
treatment has neither significant short-term benefits nor any long- 
term benefits after the drug has discontinued. It is becoming 
increasingly common to treat tumors with several drugs. In 



addition to tumor specific drugs, a generic mitotic inhibitory drug, 
which disrupts microtubules that pull the cell apart, is often used - 
since cancer cells are more sensitive to inhibition of mitosis than 
normal healthy cells. In our model, the effect of such a drug is to 
increase the death rate parameter in the equation for cancer cells. 
Further work should also include the combined effect of treatment 
of IL-27 with a mitotic inhibitory drugs. 

Methods 

Estimates of the densities of tumor cells and T cells 

Many of the parameters are based on experiments reported in 
[30]. In [30], the volume of the tumor was measured in days 
12,14,16, and 18, but the number of CD8 + T cells and 
concentrations of /10 and I y were measured only in the first 5 days. 

From Figure 5D in [30], the volume of the tumor at days 
12,14,16, and 18 were approximate 30,60,100, and 180 in mm 3 . 
Hence, 

c(12) = 30x^, c(U) = 60xA, 
c(16)=100x^, c(18)=180x^, <! 

where A is the number of tumor cells in per mm 3 . If we consider a 
simplified equation for (5) 




= hc{t)-n c c(t), 



then, for any two times t\ and ?2 5 

If we apply this formula to the 6 pairs of the numbers from (11) 
to compute k\ —fi c and take the mean value, we get 

M — n c = 0.295 /day. 

Since the half-life of melanoma tumor cells is approximate 4 
days, [47], 

\i c = In2/A / day x 0. 1 73/ day, 

and then X\ = 0.468 /day. 

In the experiments in [30], there were two kinds of tumor cells: 
J558-IL-27 which generates IL-27, and J558-Ctrl which does not 
generate IL-27. The antigen P1A onJ558 tumor cells is recognized 
by receptors TCRs on cytotoxic T cells, P1CTL. Liu et al. [30] 
used P1CTL with glycoprotein CD8 (which is called CD8 + T cells) 
to investigate the immune response for IL-27. Their P1CTL cells 
were of four different types: (i) P1CTL which can recognize J558- 
Ctrl tumor cells and generate IL-10 to inhibit tumor growth; (ii) 
IL-10" /_ P1CTL which can recognize J558-Ctrl tumor cells but 
cannot generate IL-10; (iii) P1CTL/IL-27 which can recognize 
J558-IL-27 tumor cells and generate IL-10; and (iv) IL- 
10" / "PlCTL/IL-27 which can recognize J558-IL-27 tumor cells 
but cannot generate IL-10. 

The number of tumor cells at / = 0 day (in [30]) was 5 x 10 6 
cells, and we assume (see Figure 5D in [30]) that they occupy 
volume 10mm 3 . Hence 
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c(0) = 5x 10 s cell /cm 3 



(12) 



7 27 (0) = 5 x I0 4 pg/cm 3 and I 27 (5) = 3150pg/cm\ (18) 



There is no data in [30] on the density of the tumor in day 5. 
We assume that this density is larger than c(0) but substantively 
smaller than the maximal capacity 10 9 cell / cm 3 . We take 



Taking c to be the average between the values at days 0 and 5 
(see (12) and (13)) and recalling (17), we get 



c(5) = 5.05x cell /cm 3 , (13) 
for J558-Ctrl with P1CTL or J558-IL-27 with P1CTL/IL-27, but 



c(5) = 5.1 x \$ cell /cm 3 , 



(14) 



for J558-Ctrl with IL-KTTICTL or J558-IL-27 with IL- 
10" / "PlCTL/IL-27, since the last two types of T cells do not 
generate 7io. 

From Figure 1A in [30], there were 1.1 x 10 6 P1CTL at day 1 
and 2.6 x 10 6 P1CTL at day 5; 10 6 P1CTL/IL-27 at day 1 and 
5 x 10 6 P1CTL/IL-27 at day 5. We assume that these CD8+ T 
cells occupy the volume of 10mm 3 for the first 5 days. Hence 



3750pg/cm 3 =e- 5da y x2 ^ x 5 x I0 4 pg/cm 3 



+2^x(5.025 x 10 8 )ce///cm 3 x(l- e - 5 ^ x2 ^), 



so that o<27 = 1 .5 x 10 5 pg/ cell /day 



Estimate of the parameters in (2) 

We consider a simplified version of Equation (2): 



T(l) = l.l x 10 8 ce///cw 3 ,r(5) = 2.6x \0 S cell/cm 3 , 
for P1CTL 



(15) 



— l —=s w T + a w T— ,Uir/loO)- (19) 

at hi + a io 



T(l)= 10* cell/ cm 3 , T(5) = 5 x 10* cell / cm 3 , 
forPlCTL/IL-27. 



(16) 



From [48], we have fi l0 = 16/ 'day. To estimate .?io, we consider 
the case of J558-Ctrl tumor cells, for which the term with Ixj is 
removed from (19): 



dho(t) 
dt 



--si Q T-n 10 I l0 (t). 



Estimate of the parameters in (1) 

Since IL-27 belongs to the IL-12 family, we take its diffusion 
coefficient and the degradation rate to be the same as for IL-12 
[48]: 



If T is constant, then 



7io(5) = e- 4 "io/ 10 (l) + — (1 - e- 4 "io). 

("10 



D hl = 1.25 x 10~ 3 cm 2 /day, 
fi 21 =2/ day. 



From the profile of P1CTL in Figure 3D of [30], we have 
7lo(l) = 300/>g-/c»? 3 at day 1 and 7i 0 (5) = lpg/cm 3 at day 5 and 
we take T to be the mean value of 7X1) and 7"(5) in (15). We then 
get 



In order to find 1x27, we use the simplified version of Equation 



(1): 



dhlit) 
dt 



= a.nc-n 21 In(t). 



If c is taken to be a constant, then 



727(5) = e - 5 "277 2 7(0)+ °^(1 -e- 5 "2T). 



(17) 



, / 3 -4davxl6/day , nn / 3 , S W x 1 - 8 x W S cell/cm 3 

lpg/cm = e y ' } x 300pg/cm H 



16 /day 



x(l-e 



-Adayx 



so that .vio = 8.89 x 10 i pgj 'cell /day. 

Next, we choose ffio = 5 x \Q 3 pg/cm 3 and proceed to compute 
am. We then consider J558-IL-27 tumor cells which can generate 
727. F° r simplicity, we take T to be the average between the values 
at days 1 and 5 (see (16)) and 72 7 to be the average between the 
values at days 0 and 5 (see (18): 



From Figure 1A in [30], 7 2 7(0) = 5 x \Q A pg/cm? and only 
5 — 10% of 72 7 remained in day 5. We assume that 7.5% of 727 
remained at day 5, so that 



T = 3 x 10 8 cell/cm 3 and 7 27 =2687 '5pg/ "cm 3 . 
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Then, the solution of Equation (19) satisfies 
/io(5) = g - 4 "io/ 10 (l)+— (s lo T + oc lo T T hl )(l-e- 4 "iO). 

MlO ^27 + a 10 



From the profile of P1CTL/IL-27 in Figure 3D of [30], we have 
7io(l) = 400/>g-/cm 3 at day 1 and 7io(5) = l&00pg/cm 3 at day 5, so 
that 



1800pg/W=e 



3 — Aday x 1 6 /day 



x AQQpg/cm 3 



+ —— [8.89 x 10- 8 x3x 10 s 
lb/day 

+ a m x3x 10 8 

2681 5pg/cm 3 4Ayx 16/rf a ^ 

26875^/cm 3 + 5000/)g/cw 3K ; ' 



(5_ e -4»'r) xl o 8 ce /// cw ^ 



1.3968 x 10 8 

X (l_ e -4»T) X 



(22) 



In (22), the left-hand side is close to 5 x 10 8 and the right-hand 

side is close to 5.52 x x 10 8 , while we take Ot = 

c + c T 

2x W-^pg/cm 3 and /? r = 9. 

Estimate of the parameters in (4) 

We assume that the diffusion coefficient of I y is the same as that 
of hi, namely, D y = 1.25 x W~ 3 cm 2 /day. Next we use the 
simplified version: 



dly(t) _ s y 



dt 



Sy+hl 



Hyly(t), 



(23) 



Therefore, we take aio = 1.128 x 10 4 pg / cell / day . 

Estimate of the parameters in (3) 

From [48], we have ii T = 0.3/day. For J558-Ctrl tumor cells, 
the term of ^27 in (3) drops out, and we consider a simplified 
version: 



where \i y = 2.\6/ day by [49]. For tumor cells J558-Ctrl (which 
do not generate In), (23) reduces to 



dlyjt) 
dt 



a r T-n Iy(t). 



(24) 



dT{i) c 

—y= s T— — firm, 

dt c + c T 

which, if c is constant, has the solution 

r(5) = e - 4 "7T(l)+ f T ° (l-e-^T). 
[.i T (c + c T ) 

Substituting 7X1) and 7X5) from (15), we get 

2.6 x 10* cell/ cm 3 -e-^* 03 !^ x 1.1 x \0 s cell/cm 3 (20) 



If T is constant, then the solution of (24) satisfies 



/ y (5) = e - 4 'V/ y (l)+ ^I(l- e - 4 "y). 



Taking T to be the average of 7X1) and 7X5) in (15) and taking 
Iy(i) from the curve P1CTL in the right part of Figure 3D in [30], 
we have 



3W0pg/cm 3 = e -4Avx2.i6M,>. x 250pg/cm 3 



+ o; y xl.8xl0 8 ce///W ^ _ ^4^x2.16/^ 
2.16/ day 



s T 



ji T c + c T 



x(l-e 



-4dayx0.. 



(21) 



Based on the fact that the (20) is close to 2.3 x 10 8 while the (16) 
is close to 2.3.57-^q^, we choose st = 1.3968 x 10 8 cell /cm 3 /day. 
Next, we consider 



dT(t) 



= s T - 



dt c + c T 1 + a rhi + Pt°t ho 

where the solution satisfies 



r(f) 



with Wt = 



(16), we get 



r(5) = e - 4, "7T(l) + 

l*T 

1 + a T I 2 i + P T o rho 



StC 



_(l_ e -^r) 



wt(c + ct) 
. Recalling 7X1) and 7X5) from 



so that a y = 3.72x 10 5 pg/ cell / day . We choose 
s y = 5x I0 6 pg/cm 3 . 

Estimate of the parameters in (5) 

We consider a simplified version of (5): 



dc ® 1 fA (A Il ° IA 

-— =X l c(t)- n c c(t) - r\ c c(t) - r\. 

dt l w + <j c 



Y Iy+S, 



c(t). (25) 



We choose s c = 300pg/cm 3 . In order to compute t]. f , we 
consider T cells IL-IO^^PICTL and IL-lO-^PlGTL/IL-27 
which do not generate 7io, so that 7io drops out of Equation (25): 



dc(t) 
dt 



--X l c(t)-ii c c(t)-t ly —>—c(t). 

'y i S c 
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If I y is constant, then the solution is c(t) = e 



which leads to 



r/- x l i c (0, 
f , y = [(A I -^)- 7 ln^x. 



'"/+«>' c(0) 



(26) 



Since c(5)/c(0) is close to 1 and the range of I y may vary from 
WOpg/cm 3 to 3400pg/cm 3 forJ558-Ctrl tumor cells and IL- 
10 _ / _ P1CTL T cells or from 750pg/cm 3 to I600pg/cm 3 for 
J558-IL27 tumor cells and IL-I(r/-P1CTL/IL27, we choose in 
Equation (26) 

Iy "I" Sc 



Recalling that X\ — ji c = Q. 295 / day, we take rj y = 0.6/day. 

Next, we choose a c = \50pg/cm 3 and proceed to estimate r\ c by 
considering T cells P1CTL and P1CTL/IL-27 which generate 
For (25), if I\o and I y are constants, then the solution is 

c(t) = e ( l ~*~ %1 ^~* r *™*c(0), and hence 



^ = [^1-^-^7— --In d« — f /day 



^{0.295/ day -0.6/ day x — ^— )x /l0 + ffc 



The concentration of 7 y with /jo is smaller than the 
concentration of I y where /jo is blocked [30]. We take j-fj ~ \- 
In [30], the concentration of I\o vary from \pg/cm 3 to 
2500pg/cm 3 . We take a c = \50pg/cm 3 , so that « 3.5; hence 
^ = 0.345/^. 
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